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Abstract 

A method for computing limits of quotients of real analytic functions in two variables was 
developed in [J. In this article we generalize the results obtained in that paper to the case 
of quotients q = f(x, y, z)/g(x, y, z) of polynomial functions in three variables with rational 
coefficients. The main idea consists in examining the behavior of the function q along certain 
real variety X'(q) (the discriminant variety associated to q). The original problem is then 
solved by reducing to the case of functions of two variables. The inductive step is provided 
by the key fact that any algebraic curve is birationally equivalent to a plane curve. Our 
main result is summarized in Theorem m 

In Section 4 we describe an effective method for computing such limits. We provide a 
high level description of an algorithm that generalizes the one developed in |4], now available 
in Maple as the limit/multi command. 


1 Introduction 


Algorithms for computing limits of functions in one variable are studied in [12] . Similar algo¬ 
rithms have been developed in uni and mi- Computational methods dealings with classical 
objects, like power series rings and algebraic curves, have been developed by several authors 
during the last two decades, pQ and m■ A symbolic computation algorithm for computing local 
parametrization of analytic branches and real analytic branches of a curve in ?r-dimensional space 
is presented in [2J. 

In jJ] Velez, Cadavid and Molina developed a method for analyzing the existence of limits 
lim/j- y^/o ;,) l( x t y)i where q(x, y ) is a quotient of two real analytic functions / and g , under the 
hypothesis that (a, b) is an isolated zero of g. In the case where / and g are polynomial functions 
with rational coefficients, the techniques developed in that article provide an algorithm for the 
computation of such limits, now available in Maple as the limit/multi command m- 

An alternative method for computing limits of quotients of functions in several variables has 
been recently developed in m- Their approach is completely different from ours, relaying on 
Wu’s algorithm as the main tool. 

In this article we generalize the methods presented in jJ] to the case of quotients of polynomials 
in three variables, under the same assumption that g is a function with an isolated zero at the 
point (a, b). The main idea consists in reducing the problem of determining the existence of 
limits of the form 


lim f(x,y,z)/g(x,y,z) 

{x,y,z)—^{a,b,c) 


(i) 


to the problem of determining the limit along some real variety X(q) associated to q (the discrim¬ 
inant variety of q). In order to achieve this one needs to study the topology of the irreducible 
components of the singular locus of X(q). The original problem is then solved by reducing to 
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the case of functions of two variables. The inductive step is provided by the key fact that any 
algebraic curve is birationally equivalent to a plane curve. Our main result is summarized in 
Theorem [T9] In Section 4 we provide a high level description of a potential algorithm capable of 
determining the existence of Q, and if the limit exists, it would be able to determine its value. 
Any of the Groebner Basis packages available may serve as a computational engine to implement 
such an algorithm. In Section 5 we present two examples that illustrate some the computation 
that would be needed in a typical problem of determining and computing a limit of this sort. 


2 Preliminaries 

2.1 Dimension of algebraic sets and its singular locus 

In this article we consider complex affine varieties defined by polynomials with real coefficients. 
If I is an ideal in the polynomial ring S = R[xi,..., x n ] , by X = V(I) we will denote the complex 
affine variety defined by I, i.e., the common zeros of I in C n . The dimension of X is the Krull 
dimension of the ring C®r S/I. Since S/I C C®k.S/I is a faithfully flat extension of rings, the 
dimension of X coincides with the dimension of R[xi,... ,x n \/I, the real affine ring of X. It 
is well known that if X is irreducible, defined by some prime ideal P C S, then the dimension 
of the domain R = S/P coincides with the transcendence degree of the field extension Icl 
(denoted by trdeg kL), where L denotes the fraction field of R. 

We recall the definition of the singular locus of a equidimensional affine variety. 

Definition 1. Let Y C C n be an affine variety, and let R = <C[xi,..., x n ]/I(/Y) be its ring 
of coordinates. Suppose that R is equidimensional of dimension r (i.e., ht(P) = r, for all 
the minimal primes P containing I(Y)). Let’s choose arbitrary generators fi,...,fk for I(Y). 
The singular locus of Y, denoted by Sing(Y), is the closed subvariety of Y defined by the ideal 
J = I(Y)+ the ideal of all (n — r) x (n — r) minors of the Jacobian matrix. 

Remark 2. 

1. The above criterion to determine SingfY) does not depend on the generators one chooses 
for I(Y). 

2. The singular locus SingfY) is a proper closed subvariety ofY, defined by those points p £ Y 
for which the rank of the Jacobian matrix [( dfi/dxj)(p)\ is less that n — r. 

3. dim (SingfY)) < dim(F) 

(See [7], Section 16.5 and El, Chapter I, Section 5). 

We will mainly focus in the following simple case: Suppose that X C C 3 is an affine variety of 
dimension 2 defined by a prime ideal P C R[x, y , z]. In this case P C R[x, y, z ] must be a prime 
ideal of height 1, and so it has to be principal, i.e., P = (h), where h £ R[x, y, z] is some real 
irreducible polynomial. Therefore, A' = 1 1(h). In this case Sing(A) is the complex affine variety 
defined by the ideal Is = (h,dh/dx,dh/dy,dh/dz) C R[x, y, z]. 

2.2 The discriminant variety 

The existence of lim( XiJ/ ^ f(x, y, z)/g(x, y , z) does not depend on the particular choice of 
local coordinates. Hence, after an appropriate translation we may always assume that p = (a, b , c) 


2 


is the origin, here denoted by O. Our objective is to compute 


lim J{x,y,z)/g{x,y,z), (2) 

(x,j/,z)—>(0,0,0) 

where f(x , y, z ) and g(x , y, z) are rational polynomial functions, and where g has an isolated zero 
at O. If q(x, y , z) = f{x, y, z)/g(x , y , z), we define the discriminant variety X{q) associated to q 
as the variety defined by the 2x2 minors of the matrix 


A = 


x y z 

dq/dx dq/dy dq/dz 


(3) 


Strictly speaking, the 2x2 minors of A , Xidq/dxj — Xjdq/dxi , are not necessarily polynomial 
functions. However, these minors can be written as 


Xidq/dxj — Xjdq/dxi = 
and therefore, if we let 


Xiigdf/dxj - fdg/dxj) Xj{gdf/dxi - fdg/dxi) 


fxi,xj = x i{gdf/dxj - fdg/dxj) - Xj(gdf/dxi - fdg/dxi), 

then the variety X(q) can be defined as the zeros of the ideal J = ( f x>y , f x ,z , f y ,z )• 

The following proposition states that in order to determine the existence of the limit 0 it 
suffices to analyze the behavior of the function q{x,y,z ) along the discriminant variety X(cf). 

Proposition 3. The limit lim^^^o Vi z ) exists, and equals L £ R, */ and only if for every 
e > 0 there is 5 > 0 such that for every (x,y,z) £ X(q) with 0 < \(x,y, z)\ < <5 the inequality 
| q(x, y,z) — L\ < e holds. 

Proof. The method of Lagrange multipliers applied to the function q(x, y , z) with the constraint 
x 2 + y 2 + z 2 = r 2 , r > 0 guarantees that if C r (0) = {(x, y, z) £ R 3 : x 2 + y 2 + z 2 = r 2 } then 
the extreme values of q(x, y , z) on C r { 0) are taken at those points p = (a, b, c) £ CV(0) for which 
(dq/dx(p), dq/dy(p), dq/dz(p)) = A (a, 6, c), i.e., at those points in X(q). 

Suppose that given e > 0 there is 5 > 0 such that for every (x,y, z) £ X(q) fl D* s the inequality 
| q(x,y,z) — L\ < e holds, where D* s = {(x, y,z) £ R 3 : 0 < \/x 2 + y 2 + z 2 < <5}. Let 
(. x,y,z ) £ D* s and r = \Jx 2 + y 2 + z 2 . If t(r),s(r) £ C r ( 0) are respectively the maximum and 
minimum values of q(x,y,z), subject to C r (0), then 

q(s(r)) — L < q(x, y,z)-L< q(t(r )) - L. 

As t(r) and s(r) £ X(q) fl C r ( 0) C X(q) fl D$, one sees that —e < q(s(r)) — L , and henceforth 
q{t(r)) - L < e. Thus, | q(x, y, z) - L\ < e. 

The reciprocal is obvious. □ 


2.3 Birational equivalence of curves 

We intend to reduce the problem of determining the existence of the limit 0 to a problem 
in fewer variables. In order to achieve this we will use the fact that any algebraic curve is 
birationally equivalent to a plane curve. This result follows from the following standard result: 

Proposition 4 (existence of primitive elements). Let K be a field of characteristic zero and let 
L be a finite algebraic extension of I\. Then there is z £ L such that L = K(z) (0, Page 75). 
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This immediately implies the following corollary: 

Corollary 5. Let X be an irreducible algebraic curve over a field k of characteristic zero, and 
let K be the quotient field of the ring of coordinates of X. Then for any x £ K — k which is not 
algebraic over k, K is algebraic over k{x), and there is an element y £ K such that K = k(x,y). 

The next theorem is a well known fact. Notwithstanding, we give a proof since we will need 
the explicit construction of the isomorphism denoted by p in the following theorem. 

Theorem 6. Let X be an irreducible space curve X in C 3 defined by polynomials with real 
coefficients, and such that the origin of C 3 is a point of X. Then there exists an irreducible 
affine plane curve Y C C 2 and a field isomorphism p : K(Y) —> K(X) so that X is birationally 
equivalent to Y. After removing a finite set of points Z C Y, if Yq = Y \ Z, then there is a 
morphism p : A' —> Yq such that p restricted to Xq = /i -1 (lo) is an isomorphism onto Yq. Both 
p and its inverse can be explicitly constructed. 

Proof. Suppose that X = V (P), where P C R[A', Y. Z] is a prime ideal. Since X is an irreducible 
algebraic curve dim(X) = dim(R[X, Y, Z]/P) = 1. Denote by M.(x,y,z) the fraction field of 
R[A', Y, Z]/P. Recall that dim(A) = trdeg K R(x, y , z). 

For dimensional reasons some of the variables x, y or z has to be transcendental over R. Suppose 
without loss of generality that x is transcendental over R. Corollary [5] implies that R(x) C 
M.(x,y,z) is an algebraic extension. By Proposition []] one can always find u = y + A z, for 
some A £ R(x), such that R(x,y,z) = R (x,u). Moreover, since this is true for almost all A, 
this element can be taken to be any real constant, except for finitely many choices. Define 
p : RfS 1 , T\ —> R[x, u] C M.(x,y,z) as the R-algebra homomorphism that sends S —> x and 
T —> u. Clearly p is surjective, and therefore, if J = ker(y>), there is an isomorphism of R- 
algebras p : R[S, T]/J 7* R[x, u]. Consequently, J £ R[S, T] is a prime ideal. Denote V(J) by 
Y. The last isomorphism induces a field isomorphism p : R(Y) = R(x, u) —> R(x, y, z ) defined as 
p(x) = x , p{u) = y + A z. Therefore, dim(F) = dirn(A) = 1. Hence, Y = V(J) is an irreducible 
algebraic plane curve which is birationally equivalent to X. 

The morphism p : R(F) —> R(A) induces a morphism of varieties p : X —» Y given by p(a, b , c) = 
(a, b + Ac). 

Notice that since (0,0,0) £ A, then, obviously, (0,0) £ Y. Since Y is an irreducible plane 
curve, J must be a height one prime ideal. Thus, J = (h), for some h(X, U ) £ R[X, U\. 

We can assume that the polynomial h(a, U ) obtained by replacing the variable A by a £ C is 
not identically zero: If h(a , 17) = 0 we would have h{ A, U) = (A — a) m t( A, U ), with t(a, U) ft 0. 
But (0,0) £ Y implies a = 0, and henceforth h(X,U) = X m t(X,U). Thus, h(X,U) = A or 
h(X,U) = t(X,U), since Y is irreducible. Finally, we note that h(X,U) = X contradicts the 
fact that x is transcendental over R. 

On the other hand, since x is transcendental over R, by Corollary [5] the extension R(x) C 
R(x)(it) is algebraic. Therefore, since y £ R(x, u) one can write y as: 

a 0 (x) a i(x) , a r (x) r 

y — , 7 \ + 1 / \ u ”1-+ 7 / \ u ; 

b 0 {x) bftx) b r {x) 

where r is smaller than the degree of the field extension [R(x)(u) : R(x)]. Taking b(x) = 
bo(x) ■ ■ ■ b r {x) we can rewrite the last equation as 

c 0 (x) + c 1 (x)u-\ - + c r (x)u r . 

v =- W) -’ (4) 

for certain cftx). Therefore, we have y = fi(x,u)/gi(x) and z = f 2 {x,u)/g 2 (x). Consider 
Z = {(a, b) € Y : gfta) = 0 or ( 72 (a) = 0}, which is a Zariski closed subset of Y. 
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Let us see that Z is a finite set. Indeed, the polynomials g\ and g 2 have finitely many roots. 
Therefore, if 01 ,..., a*, £ C are these roots, for each at, (a*, b ) £ Y if and only if h(a,, 6) = 0, 
where Y = V(J) with J = ( h ). Notice that the polynomial t(U) = h(a,i, U ) £ C[C/] has finitely 
many roots. Hence, there are only finitely many elements (a,, 6) with gi{af) = 0 or g 2 (ai) = 0, 
and such that h(at, b ) = 0. Thus, we conclude that Z is finite. 

Consider the open subset Y 0 = Y \ Z of Y. Let Xo = /r _1 (Y 0 )- Define r : Yo ► X 0 as r(d, e) = 
(■ d , ; ^g^d) ) • This last morphism induces an R-algebra homomorphism ip : R(X) —>• Oy(Yo) 

given by ip(x) = s, ^(y) = fi(s,t)/gi(s), and ip(z) = f 2 (s,t)/g 2 (s). Clearly, 

/ / \ fl(x,u) f 2 {x,u) 

tpo ip{x)=x, <p o ipyy) = - -—=y, tpo ip{z) = - -— = z. (5) 

9 iW 52(a:) 

Therefore, p o ip = Id^x), and consequently <p o ip\x 0 ■ Ox\x 0 —> Oy\y 0 is th e identity. On 
the other hand, ip o <p(s) = ip{x) = s and ip o ip(t) = ip{u). By © we have <p o ip(u) = u and 
tp(t) = u, which implies that t = ip(u), since p is injective. Hence, ip o ip(t) = t and therefore 
ip o p\ Yo : Oy\y 0 —> Ox\x 0 is the identity. Hence, ip : Ox\x 0 —> Oy\y 0 is the inverse of the 
morphism <p : Oy\y 0 —> Ox\x 0 - Thus, the homomorphism r : lb —i> Xo induced by ip is the 
inverse of p : X 0 —> Lo- 

Finally, it is clear that the morphism p : Xq —> Yq sends the real part of Xo into the real part of 
lo, and since p^ 1 = t : Yq —> Xq is determined by the polynomials /i,/ 2 ,<?i and g 2 , which are 
all real polynomials, then p~ 1 = r also sends the real part of Yq into the real part of A'o- □ 

Remark 7. Xo is obtained from X by removing finitely many points. 

Proof. In fact, a point (a, b, c) £ X' does not belong to X 0 iff p(a, b , c) = (a, b + Ac) ^ lo, he., iff 
(a, 6+Ac) £ Z. But Z is finite, and therefore there are only finitely many choices for a and (6+Ac) 
such that (a, b + Ac) ^ Ifo Fix any values for a and for r]=b + Ac. If fi(x , y, z),..., fk(x,y, z) are 
generators for P then, clearly, /)(a, — Ac, c) = 0. But each polynomial gt(z) = fi{a, p — A z, z) 

can only have finitely many roots. This proves the claim. □ 

This Remark tells us that the problem of determining (and computing) the limit of a function 
along the varieties X and Y is equivalent to the same problem when one approaches the origin 
along Xo and To- 

2.4 Groebner bases 

In this section we collect some basic properties and results on Groebner bases and Elimination 
Theory that will be needed later for the development of an algorithm that computes ([]]). The 
main reference for this section is [7], Chapter 15. 

By S = K[x\,..., x n ] we denote the polynomial ring in n -variables with coefficients in a field 
I\. We denote the set of monomials of S by M. By a term in S is meant a polynomial of the 
form cm, where c yf 0 £ K and m £ M. 

Definition 8. A monomial order in S is a total order on M satisfying nm\ > nm 2 > m 2l for 
every monomial n ^ 1, and for any pair of monomials m\ and m 2 satisfying m\ > m 2 . 

Every monomial order is Artinian which means that every subset of M has a least element. 
For a fixed monomial order > in S', the initial term of p £ S is the term of p whose monomial 
is the greatest with respect to >. It is usually denoted by in(p). Given an ideal I C S, its ideal 
of initial terms, in(J), is defined as the ideal generated by the set {in(p) : p £ /}. 
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Definition 9. Let I C S be any ideal, and fix a monomial order in S. We say that a set of 
elements {/i,..., fk} of I is a Groebner basis for I iff in(I ) = (m(/i),..., in(fk))- 

We list some basic facts about Groebner bases. 

Remark 10. 1. The set of monomials not in the ideal in(I) forms a basis for the K-vector 

space S/I. 

2. There always exists a Groebner basis for an ideal I C S. As S is a Noetherian ring, 
the ideal I is finitely generated, let’s say, I = (/i, ■■■,/&)■ Consider the ideal J = 
(m(/i),..., in(fk))- If J = in(I) then {/i,..., fk} is a Groebner basis for I. 

3. If {/i,..., fk} is a Groebner basis for I then I = (/i, ...,/*,)■ 

f. There is a criterion that allows to compute algorithmically a Groebner basis for an ideal 
I C S. This criterion is known as Buchberger’s algorithm m Page 332). 

5. Let I, J be ideals of S such that I C J. If in(I) = in{J) then I = J. 

An example of a monomial order is the lexicographic order, defined in the following way: Fix 
any total order for the variables, for instance x\ > > • • • > x n , and define x'fixff ■ ■ ■ xf n > 

x b fix b fi ■ ■ ■ x if for the first j with aj 7 ^ bj one has aj > bj. (The lexicographic order will be the 

monomial order that we will use in this article.) 

Now we discuss a basic result that will be needed in Sections 4 and 5. 

Let / be an ideal of the polynomial ring K[x 1 ,... x n , yi,..., y s \. Given a Groebner basis for 
I we want to compute a Groebner basis for I D K[x 1 ,..., x n ]. For this, we have to introduce the 
notion of an elimination order: 

Definition 11 . A monomial order in K[x\,... x n , y \,..., y s \ is called an elimination order if 
the following condition holds: f £ K[x \,... x n , yi ,..., y s ] with in{f) £ K[x\,...,x n ] implies 
f £ K[. Ei, .. .,x n ). 

Lemma 12. Let I C K[x 1 ,... x n , y \,..., y s \ be an ideal, and let B = {/ 1 ,..., fk} be a Groebner 
basis for I with respect to an elimination order. Assume that fi,...,ft with t < k are all 
elements of B such that fi,..., ft £ K[x 1 ,..., x n \. Then {/ 1 , • • ■, /t} is a Groebner basis for 
i n K[x 1 ,.. .,x n ]. 

Proof. (See J7], Page 380). □ 

Remark 13. Suppose that (p : K[x 1 ,..., x n ] K[yi, ■ ■ ■, y s ]/J is a ring homomorphism defined 
as tp(xi) = fi. Consider Fi £ K[yi ,... ,y s ] such that Fi = fi in K[yi,... ,y s \/J, and define the 
ideal I = JT + ( F\ — X\,... ,F n — x n ) C T, where T = K[x 1 ,... ,x n ,y\,... ,y s \. Then ker ip = 
IC\K[x 1 ,..., x n ]. Therefore the above lemma implies that ker p can be computed algorithmically. 

Proof. (See [7|, Page 358). □ 

3 Reduction to the case of functions of two variables 

Let q(x,y,z) = f(x,y,z)/g{x,y,z) be the quotient of two polynomials. We recall ('Section [ 2 . 21 ) 
that the discriminant variety associated to q, X(q) C C , is the affine variety defined by the 
2x2 minors of the matrix we denoted by A. As a variety, X(q) may be decomposed into its 
irreducible components in C 3 , let’s say X(q) = X\ U X 2 U • • • U X n . 

We are only interested in those components that contain the origin. These will be called the 
relevant components. Suppose these are Xi, X 2 ,..., Xk, k < n. We consider three possible cases: 
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1. dim Xt = 0: In this case, if Xj = V{Pf), then R[a :,y,z\/Pi is a field and Xj is just the 
origin {O}. Hence, X, does not contribute to any trajectory in R 3 that approaches O, and 
can be discarded. 

2. dim X,; = 1 : In this case Xj is an irreducible algebraic curve. 

3. dim Xj = 2: In this case Xj is an hypersurface, i.e., Xj = V(Pi ), where Pi is a principal 
ideal. 

We only have to study Cases 2 and 3. 

We deal first with the case of an irreducible space curve in C 3 . Let us see that the problem 
of determining the limit of q(x,y,z ) along X', as well as its computation can be reduced to the 
case of a real plane curve, a question already addressed in Jj . 

By Theorem [Gl there is a plane curve Y which is birationally equivalent to X, and therefore 
a local isomorphism y : Xo —> lo, where Xo and Id are as in Theorem [Gj There we observed 
that the existence of the limit of q(x, y, z) as (x, y, z) —> (0, 0,0) along X 0 is equivalent to the 
existence of the limit of q o y~ x as (it, v) —> (0,0) along Y 0 . Thus, 

lim q(x,y,z) = lim qo y~ l {u,v). ( 6 ) 

{x,y,z)->0 (u,v)->0 

(x,y,z)£X 0 ( u,v)gY 0 


Summarizing: 

Proposition 14. Let X C C 3 be an irreducible component of X(q) of dimension 1 containing 
O. Let y : Xo —> Yq be the local isomorphism defined in Theorem [ 6 J Then, the limit ofq(x, y,z ) 
as (x, y,z) O along X exists if and only if exists along the irreducible plane curve Y as 
(u,v) —> (0,0). The corresponding limits are related by ffj). 

In the sequel we will denote by Xj and Y, the two open subsets X ' 0 C X and To C Y defined 
in Theorem [ 6 ] For the purpose of analyzing the limit along the space curve X' it is only necessary 
to consider those cases where the real trace of the birationally isomorphic curve Y turns out to be 
a plane curve containing the origin. By yx t we denote the corresponding isomorphism between 
Xj and Y.j already constructed. 

Now we analyze Case 3. This is a lot more subtle, and requires a careful analysis of the 
topology of the corresponding two dimensional component. A key ingredient is a celebrated 
theorem of Whitney m about the number of connected components of an affine algebraic variety. 
In the following discussion we will show how one can reduce the analysis of the 2-dimensional 
irreducible components to Case 2. 

Suppose that we have a rational function q(x,y,z ) = f(x,y,z)/g(x,y,z) defined on an irre¬ 
ducible hypersurface X = V(h), where h is a real polynomial function of three variables and q 
has an isolated zero at 0. Let S = Sing(X) be the singular locus of X. By Remark [2] S must 
be a variety of dimension strictly less than two. Hence, if S contains the origin, the limit of q as 
(x, y,z) —» O along S can be computed as in Case 2. 

Now, we restrict our analysis to the nonsingular locus of X, that we denote by A f = X \ S. 
Without loss of generality we may assume that JV contains the origin, otherwise all of their 
components would be irrelevant. 

Assume O £ AT, and define a family of real ellipsoids E r = {(x,y,z) € R 3 : Ax 2 + By 2 + 
Cz 2 — r 2 = 0}, A, B, C > 0, r ^ 0. By p r (x,y,z) we will denote the quadratic polynomial 
Ax 2 + By 2 + Cz 2 - r 2 . 
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Definition 15. Let X = V(h) C C 3 and E r = {(x, y, z) G R 3 : Ax 2 +By 2 +Cz 2 —r 2 = 0}, r ^ 0 
as above. The critical set C r (q) will be the set of all real points in E r C\X where q(x , y , z) attains 
its maxima and minima. The union U r >o C r (q) of all critical sets will be denote by Critx(g)- 

Since each E r fl X is a compact set, and by hypothesis O is an isolated zero of q, the set 
Critx(g) is a well defined subset of X. 

We need the following analogue of Proposition [3] 

Proposition 16. The limit lim^^^^o q{x , y, z) along X exists and equals L if and only if for 
every e > 0 there is 5 > 0 such that for every 0 < r < 5 the inequality \q(x, y,z) — L\ < e holds 
for all (x, y, z) € C r . 

Proof. The proof follows identical lines as in Proposition [3] One just have to notice that each 
point in the critical set must lie in some E r , since p = (a.b.c) is obviously contained in E r , with 
r = V Aa 2 + Bb 2 + Cc 2 . " □ 

Our objective is to determine Crit.v(g)- We can decompose this set as the union of CritAf(g) = 
Critx(g)nA/’ and Critx(g)n«S. Since Critx(g)05 C S , and the limit along S can be determined 
as in Case 2, we just have to focus on Crit jq(q). 

First, we want to determine the nonsingular part of Crit/v(g) by using the method of Lagrange 
Multipliers, as in 0. For this we define X = P(3) C X to be the zero set of the ideal 5 generated 
by h and the determinant: 


d{x,y,z) = 


dp r /dx dp r /dy dp r /dz 
dh/dx dh/dy dh/dz 
dq/dx dq/dy dq/dz 


As the points of X already satisfy X7q(x,y,z) = X(x,y,z) (where Vg denotes the gradient of 
g), and since Xp r (x,y, z) = (2Ax,2By,2Cz), the affine variety X must be defined by the ideal 
generated by h and by the determinant: 


D{x,y,z ) = 


Ax By Cz 

x y z 

dh/dx dh/dy dh/dz 


That is, X = V(D,h). This variety is precisely the set of regular points of X that are critical 
points of q. 

Proposition 17. (Notation as above) Let us assume O € Af. Then it is possible to choose (in a 
generic way) suitable positive constants A, B and C such that the height of the ideal 3 = (D, h) 
in the polynomial ring C[x,y,z] is greater than one, and consequently dimX < 2. 

Proof. It suffices to show that for a suitable choice of positive constants A, B, C there is at least 
one point p ^ O in M such that D(p) ^ 0. 

First, let us see that there is at least one point p £ Af different from the origin such that 
the gradient of h does not point in the direction of p, i.e., such that Vh(p) ^ A p, for all A G R. 
Indeed, suppose on the contrary that for every p £ Af there existed A (p) ^ 0 such that S7h(p) = 
A (p)p. Since each p is a regular point of X , one must have Xh(p) ^ 0. Hence, after making 
an appropriated change of coordinates that fixes O (a rotation, and then a homothety) we may 
assume without loss of generality that dh/dz{ 0, 0, 1) ^ 0, and that p = (0, 0,1). By the implicit 
function theorem there would exist Uq C R 2 , a neighborhood of (0,0), and a smooth function 
u(x,y) in Uq such that u(0, 0) = 1, and h(x,y,u(x,y)) = 0, for all (x,y) £ Uq. Since A7h(p) = 
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A (p)p, one must have dh/d x( 0 , 0 , 1 ) = dh/dy(0,0, 1 ) = 0 , and consequently du/dx(0,0) = 0 = 
du/dy{ 0 , 0 ). 

Let W p be the graph W p = {(a;, y, u(x, y)) : (x, y) € Uq}. For any t S W p , the normal vector 
at t is given by 

( 'U‘x j 'Uyi 1 ) 


»(*) = 


ul+ul + 1 


Henceforth, if y(t) = A(t)/||V/i(t)|| one has that Vh( t) = /x(t)||V/i(t)||t, and consequently n(t) 
can be written as 

(x,y,u(x,y)) 

n(t) = 


\Jx 1 + y 2 + u 2 (x, y) 


From this, we deduce: 


1 

yjul+ul + 1 
~u x 

yjul+ul + 1 


u{x, y) 

\/x 2 + y 2 + u 2 (x, y) 
x 

\Jx 2 + y 2 + u 2 (x,y) ’ 


and 

~u y = _ y _ 

^Jul + zi 2 + 1 V x2 + V 2 + u 2 (x,y) 

This implies u x = —x/u(x,y), and u y = —y/u(x,y). Hence, u(x,y) = \]l — x 2 — y 2 1 since 
u(0,0) = 1. We conclude that W p would be a neighborhood of p in JV which is part of a sphere 
centered at the origin. But on the other hand, a theorem of Whitney asserts that M can only 
have finitely many connected components (see [20]). Then this would imply that M could not 
contain the origin, a contradiction with our assumption. 

Therefore, we may assume there exists a point p ^ O in AT such that V/i(p) ^ A p, for all 
A ^ 0. After applying a rotation (if necessary) we may also assume that a, b , c are all nonzero. 

After those preliminaries it becomes clear how to choose positive constants A, B and C such 
that the determinant 


Aa Bb Cc 

a b c 

dh/dx(a,b,c) dh/dy(a,b,c) dh/dz(a,b,c) 

does not vanish: The vectors V/i(p) and p = ( a,b,c ) generate a plane H , since they are not 
parallel. Therefore, it suffices to choose any point (a,/?, 7 ) outside H and such that A = a/a, 
B = fj/b, and C = c /7 are positive. □ 

As before, for the limit lim ^ iViZ )^o dix, y, z ) to exist along X it is necessary that it exists 
along any real curve that contains O. In particular, the limit along each component of X must 
exist, and all theses limits must be equal. By Proposition 17, dim(X) < 2, and henceforth we 
can reduce this last question to cases 1 and 2 . 

Let 3 be the affine variety defined by the ideal generated by h and by the minors 2 x 2 of the 
matrix 

Ax By Cz 
dh/dx dh/dy dh/dz 
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The set 3 H E r (~1 AT defines the locus of those real points where E r and A f do not intersect 
transversely. Outside this set, E r D Af is a 1-dimensional manifold (see [15] . Page 30) that we 
shall denote by E. Clearly, the vanishing of these two by two minors forces the vanishing of the 
determinant D(x,y,z). Henceforth, 3 C X, and consequently dim(3) < 2, by Proposition H71 
Again, for the existence of the limit lim ( Xi y iZ )^o q( x ,y, z) it is required, in particular, its exis¬ 
tence along any relevant component of 3, and consequently the problem reduces again to cases 
1 and 2. This takes care of the subset of Critv(g) inside 3- 

As for those points in Critv(g) that lie outside 3, we notice that they are contained in the 
1-dimensional manifold E. Then they must be part of X, since this variety is precisely those 
regular points where q attains an extreme value. Thus, the points in Crit/v(g) that lie outside 3 
must be contained in X. Once again, we have reduced the problem to cases 1 and 2. 

The following proposition summarizes this discussion: 

Proposition 18. Let X be a relevant irreducible component of dimension 2 of the discriminant 
variety X(q). Consider S, X, and 3 as defined above. Then, the limit of q(x,y, z) as ( x,y,z) —> O 
along X exists, and equals L, if and only if, the limit of q(x, y , z) as (x, y,z) —► O exists and 
equals L along each one of the components of the curves S, X, and 3- 

We are now ready to state our main result. 

Theorem 19. Let q{x,y,z) = f(x,y,z)/g(x,y,z), where f and g are rational polynomial func¬ 
tions, and where g has an isolated zero at the origin. Let X(q) be the discriminant variety 
associated to q. Denote by {X±,...,X^} the relevant irreducible components of dimension one 
of X(q), and by {Xfc + i,... ,X n } the relevant irreducible components of dimension two of X(q). 
Then, the limit of q as (x, y,z) —> O exists, and equals L, if and only if the limit of q(x, y, z) as 
(x, y, z) —> O along X j exists, and equals L, for all i = 1,2,..., n. Moreover: 

1. For the components Xi, i = 1,2, ...,k, the limit of q(x,y,z) as ( x,y,z ) —> (0,0,0) along 
Xi is determined as in Proposition \14\ 

2. For the components Xj, j = k + 1,..., n, the limit of q(x, y, z) as (x, y, z) —> (0,0, 0) along 
Xj is determined as in Proposition 1 1 #1 


4 A high level description of an algorithm for computing 
the limit 

Let q{x,y,z) = f(x,y,z)/g(x,y,z), where / and g are polynomial functions of three variables 
with rational coefficients, and g has an isolated zero at the origin. Consider X(q), the discrimi¬ 
nant variety associated to q. We have to decompose X(q) into irreducible components, and then 
choose only those irreducible components {Xi, ..., X n } that are relevant. 

The algorithm has to deal with two different cases: 

• Dl: The component Xi has dimension 1. Then as observed before, Xj is birationally 
equivalent to an irreducible plane curve Y). Let us denote by C(x ,y,z) the fraction field 
of the ring of coordinates C[X,Y, Z\/I(Xj) of Xj. As we already noticed we may always 
assume that x, y, z are transcendental elements over C: If, for instance, x were algebraic over 
C, then there would exist a polynomial P(X) £ C[A] such that P(x) = 0. This is equivalent 
to saying that P(X) £ I(Xi). Suppose we write P(X) = (X — a\){X — a 2 ) • • • (X — a n ) 
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in C[X]. Since J(Xj) is a prime ideal, some linear factor X — aj must belong to /(X,;). 
But as X,; contains the origin, we must have ay = 0. Hence, we could write I(Xi ) = 
(X, hi(Y, Z ),..., h m (Y, Z)), where hk(Y, Z) £ C[Y, Z], for k = 1,2,..., m. If we denote 
by I the ideal {h\{Y, Z ),..., h m (Y, Z)), and by X- = V(I) C C 2 the affine variety defined 
by I, then the limit of q(X, Y, Z) as (X, Y, Z) —► O along Xi is the same as the limit 
of q(0,Y,Z) as (Y,Z) —> (0,0) along X'. But the existence of this limit, as well as its 
value, can be computed using the algorithmic method developed in [3j. By Proposition 
Eland Corollary [5] we know that if x is transcendental over C there exists A £ C(x) such 
that C(x,y,z) = C (x,u), where u = y + A z. Also, by Theorem [Gj if we consider ip : 
C[X, U ] —> C[X, Y, Z\/I(X) defined by ip(X) = x, ip(U) = y + A z, then kerip defines the 
irreducible plane curve Y that is birationally equivalent to X. As we observed in Section 
12.41 ker<^ = (I(W)T + (U — (Y + A Z))) D C[X, U], where T = C[X, U, Y, Z] is computable. 
On the other hand, the ring homomorphism C[X, U]/ kery> C [X,Y, Z]/I(X) induces 
an isomorphism of fields ip : K(X) —> K(Y). As we showed in the proof of Theorem El 
since y,z £ C (x,u) then one must have y = fi(x,u)/gi(x), and 2 = / 2 (x, u)/g 2 {x), for 
some f t , f '2 , gi , and g 2 with real coefficients. In that same proof we noticed that the local 
isomorphism /i : Xq —> Yq is determined by those polynomials. 

By Proposition fldl computing the limit of q(X,Y,Z) as (X,Y,Z) —»• O along X- L is equiv¬ 
alent to computing the limit of q o p^}(X,U), as (X,U) —> (0,0) along Y), and this last 
limit can be dealt with using the algorithm developed in [3j. 

• D2: Suppose that dimXj = 2. Then X,; is an affine variety defined by a principal ideal 
I{Xi) = ( h ). For random positive values A , B and C the algorithm computes the height 
of the ideal J = ( D,h ), where 


D = 


Ax By Cz 

x y z 

dh/dx dh/dy dh/dz 


As we saw in the reduction to plane curves, there always exist positive constants A , B and 
C such that ht(3) > 2. Since dim(X) < 1, then X = V(3). Henceforth, one can compute 
the limit of q{x , y , z) as (x, y,z) O along X using the prescription in Dl. 

Since S = Sing(X), the affine variety defined by the ideal (h, ^)) must be a proper 

subset of X. Then S is also an algebraic curve, and once again we can compute the limit 
of q(x, y , z) as (x, y,z) —» O along S using Dl. 

Now, the affine variety 3 defined by the ideal generated by the minors 2 x 2 of the matrix 

Ax By Cz 

dh dh dh 

. dx dy dz . 

and the polynomial h, has also dimension less than 2. Hence, the limit of q(x,y,z ) as 
(. x,y,z ) —> O along 3 is also computed using Dl. 

• Finally, if the limit of q(x, y 7 z) as (x, y, z) —»• O along each relevant irreducible component 
of X(q) of dimension one exists, and equals L , one says that the limit of q(x,y,z ) as 
(, x,y,z ) — > O is L. Otherwise, one says that this limit does not exists. 
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5 Examples 

5.1 Example 1 

Suppose that we want to compute the limit: 


lim YX ~ ZY + ZX 
(x,y,z)™( o,o,o) X 2 + Y 2 + Z 2 

Let q(X, Y, Z) = YX — ZY + XZ/X 2 + Y 2 + Z 2 . We illustrate the necessary computations, 
carried out in the program Maple. 


1. Using the command PrimeDecomposition(X(q)) one gets the irreducible components of 
X(q): 

V({Y-X + Z)), V{{X 2 +Y 2 + Z 2 )), V((X + Y,Z- 2X)), 

V{{X + Y,Z + X)) and V((X + Y, Z 2 + 2X 2 )). 


By using the command HilbertDimension(Q) one can see that the irreducible component 
V(X + Y, Z + X) has dimension 1. 

Let us see that for A = 1, C (x,u) = C(x,y,z), where u = y + z. Here, C(x,y,z) denotes 
the fraction field of the ring of coordinates of the variety V(X + Y, Z + X). Consider the 
ideal I = (X + Y, Z + X)T + (U - (Y + Z)), where T = C[A, Y, Z, U]. The command 
Eliminationldeal (I ,{U,X}) generates the ideal J = (2A + U). On the other hand, the 
command Basis (I ,plex(Z, Y,U) ) gives us a Groebner basis for I respect to the lexico¬ 
graphic monomial order, with Z > Y > U. In this particular case we obtain the following 
basis: {2X + U, Y + X, Z + X}. From this basis we deduce that y = — x and z = —x are ele¬ 
ments of C(x, u). Therefore, C(x, u) = C(x, y, z), and consequently the ideal J = (2A + U ) 
defines an irreducible plane curve which is birationally equivalent to V(X + Y, Z + X). 
Also, y = —x and z = —x determine the isomorphism p : V(2X + [/)—>■ V(X + Y, Z + X). 
Therefore, the limit of q(X, Y, Z), as (X, Y, Z) (0, 0,0) along V(X + Y,Z + X), is equiv¬ 
alent to the limit of q o p(X, U ) as (A, U) —> (0, 0) along V(2X + U). This latter limit can 
be computed using the algorithm developed in [5]. However, in this case it is easy to see 
directly that the value of the limit is —1/3, since q o p(X, U) = q( A,—A,—A) = —1/3. 
Therefore, the limit of q( A, Y, Z) as (A", Y, Z) —> (0, 0,0) along V(X + Y, Z + X) is —1/3. 
Let h(X,Y,Z) = Y — X + Z. One may choose A = 1, B = 2, and (7 = 1. Using the 
command HilbertDimension(P) , with P = ( f,h ), where 


X 

2 Y 

z 

X 

Y 

z 

dh 

dh 

dh 

dx 

8Y 

dZ 


one obtains that the variety defined by the ideal P = (/, h) has dimension 1. In this case 
V(P) = V(—XY — YZ,Y ~ X + Z). Using again the command PrimeDecomposition(P) 
one obtains the irreducible components of the variety V(P): V(P) = V(Y,—X + Z) U 
V(2X — Y. Y + 2Z), where each of these components has dimension 1. Therefore, one 
just needs to compute a limit along irreducible algebraic curves (again, using the main 
algorithm of [4]). For the variety V(Y + A, Z + X) we may follow an analogous procedure. 
It is not difficult to see that the limit of q( A, Y, Z), as (A, Y , Z) — > (0, 0,0) along the variety 
V(Y, — A + Z) is equal to 1/2. 

Hence, we conclude that 


lim 

(X,Y,Z)->( 0 , 0 , 0 ) 


YX - ZY + ZX 
X 2 +Y 2 + Z 2 
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does not exist. 


5.2 Example 2 

1. We want to compute the limit: 


X 2 YZ 

lim - 

(x,y,z)->(o,o,o) X 2 + F 2 + F 2 

Let q(X, F, Z) = X 2 YZ/X 2 + Y 2 + Z 2 . 

Using the command PrimeDecomposition(X(q) one obtains the irreducible components 

of X(q): 

2. The irreducible components of dimension 1 are: Xi = V(X,Y), X 2 = V(X,Z), X 3 = 
V(Y,Z), X 4 = V(X,Z-Y), X 5 =V(X,Z + Y), X 6 = V(X,Y 2 + Z 2 ), X 7 = V(X,3Z 2 - 
Y 2 ), X 8 = V(X,3Z 2 +Y 2 ), X 9 = V(Z,X 2 + Y 2 ), X 10 = V(Z + Y,-2Z 2 + X 2 ), X n = 
V(-Z + Y,-2Z 2 + X 2 ), and X 12 = V(-2Z 2 + X 2 ,3Z 2 + Y 2 ). 

3. The irreducible components of dimension 2 are: X 13 = V(X), and X 74 = V(X 2 +Y 2 +Z 2 ). 

4. There is an irreducible component of dimension 0: X 15 = V(X, Y, Z). 

We know that each irreducible component of dimension 1 is birationally equivalent to 
an irreducible plane curve. Now, if any of the variables X, Y or Z appears in the ideal 
that defines the corresponding irreducible component, then one can easily see that such 
component in C 3 is actually contained in C 2 . Henceforth, it is already a plane curve, and 
one can use the main algorithm of [3] to compute these (two variable) limits. Hence, one sees 
that the limits along the varieties X\ = V(X,Y), X 2 = V(X,Z), X 3 = V(Y,Z), X 4 = 
V(X,Z-Y), X 5 = V(X,Z + Y), X 6 = V(X,Y 2 + Z 2 ), X r = V(X,3Z 2 -Y 2 ), X 8 = 
V(X, 3 Z 2 + Y 2 ), and Xg = V(Z, X 2 + Y 2 ) are equal to zero. 

Now we discuss the limit along the other irreducible components of dimension 1. 

5. For X 10 = V ( Z+Y , —2Z 2 +X 2 ): We noticed in the proof of the Primitive Element Theorem 
that for almost all A £ K, M.(x,y,z) = M(x, u), where u = y + A z. In this case, one could 
take A = 2. Let I = (Z + Y, — 2Z 2 + X 2 , U — (Y + 2Z)) C M[X, Y, Z, U], with the command 
EliminationIdeal(I, X, U) one gets the plane curve V(2U 2 — X 2 ), which is birationally 
equivalent to Xio- On the other hand, by using the command Basis (I,plex(Z,Y,U)) one 
computes the basis {2C7 2 — X 2 , Y + U, —U + Z}. From this basis we deduce that y = —u 
and z = u as elements of R(x, u). Thus, the limit along the component X 10 is the same as 
the limit of q(X , — U, U ) along the irreducible plane curve V(2U 2 — X 2 ). This latter limit 
can be calculated using the the main algorithm of . In this case we obtain the value zero. 

6 . For Xn = V(—Z + Y,—2Z 2 + X 2 ), with A = 1 and following the same procedure, 
i.e., defining the ideal I = (-Z + Y,—2Z 2 + X 2 ,U — (Y + Z)) and then computing 
EliminationIdeal(J, X, U ) and Basis (I ,plex(Z, Y, [/)), one obtains the irreducible plane 
curve V(U 2 — 2X 2 ), which is birationally equivalent to Xu, as well as the basis {I/ 2 — 
2X 2 , — U + 2Y, —U + 2Z}. From this basis one deduces that y = u/2 and 2 = u/2, and 
therefore the limit along the component Xu is the same as the limit of q(X , U / 2 , U/2) 
along the plane curve V(U 2 — 2X 2 ) which is again a limit in two variables, and can be 
computed using the methods of [4] . In this case the limit is also zero. 
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7. For X 12 = V(—2Z 2 + X 2 ,3Z 2 + Y 2 ), with A = 1, by defining the ideal 

I = ( -2Z 2 + X 2 , 3Z 2 + Y 2 ,U - (Y + Z)), 

and then computing EliminationIdeal(/, X, U) and Basis (I ,plex(Z,Y,U)), one obtains 
the irreducible plane curve V{4X 4 + 2X 2 U 2 + U 4 ), which is birationally equivalent to X 12 , 
as well as the basis 

{4X 4 + 2 X 2 U 2 + U 4 , - U 3 - 4UX 2 + 4X 2 Y, 4ZX 2 + U 3 }. 

From this we deduce y = u3 \ 4 2 X , and z = and therefore the limit along the 

component X 12 is the same as the limit of q(X, u3 ~j[ffl x , along the plane curve 
V(4X 4 + 2X 2 U 2 + U 4 ), which is again a limit in two variables. This limit is also zero. 

8 . Now, the components of dimension 2 are V(X 2 + Y 2 + Z 2 ) and V(X). The first component 
is precisely the set of points where the rational function q is not defined, and consequently 
can be discarded. Since the variable X appears in the ideal defining the second variety the 
limit clearly must be zero. 

We conclude that 

V 2 VZ 

(x,y,z)™(o,o,o) X 2 + Y 2 + Z 2 ! 

since it is zero along each of the irreducible components of the discriminant variety. 
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